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CN ■ Abstract 
^ , 

CN i Within the framework of finite temperature field theory this paper 

' discusses the shear viscosity of hot QED plasma through Kubo formula 

■ at one-loop skeleton diagram level with a finite chemical potential. 

^ i The effective widths (damping rates) are introduced to regulate the 

Q^' pinch singularities. The finite chemical potential, which enhances the 



Oh! 



contributions to the shear viscosity from the electrons while suppresses 



D , those from the photons, finally gives a positive contribution compared 



X 

- . 1 Introduction 



to the pure temperature environment. The result agrees with that 
from the kinetics theory qualitatively. 



Transport properties of relativistic plasma are of great interest both exper- 
imentally and theoretically. Taking the signal of quark gluon plasm(QGP) 
for example, the so-called 'strong-coupled matter' formed in the relativis- 
tic heavy ion collider(RHIC) was well described by an ideal hydrodynamics 
model in the region of pt < 2GeV to fit the data of the elliptic flow V2 [1,2]. 
But it is pointed out that the over-prediction in high-pt region of such model is 
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due to the neglect of dissipation and viscosity [3,4]. Thereby the transport co- 
efficients need to be considered. In addition, the non-ideal fluid properties of 
viscosity and dissipation may influence some other important equations and 
quantities such as the equation of state(EoS), the formation (thermalization) 
time and the sound attenuation length etc. [4-7]. 

From the theoretical points of view, two formalisms are usually developed 
to calculate the transport coefficients. One is kinetic theory which starts 
from expanding the distribution function near the local equilibrium when 
dealing with the Boltzmann equation. The transport coefficients, especially 
the shear viscosity, were discussed by many authors in this formalism [8-14]. 
The other way of obtaining the transport coefficients is to implement the 
Kubo formulae within the framework of the finite temperature field theory. 
The Kubo formulae, which work in both weak coupling and strong coupling 
regime, have been applied for computing the transport coefficients of strong- 
coupled theory in lattice Monte Carlo simulations [15-17], and for analytical 
calculation in scalar theory [18-20] and gauge theory [21-24,27,28], early 
with relaxation time approximation and later with ladder diagram resum- 
mation and large Nf expansion technique. However these two approaches 
are not irrelevant. Their inner relations are attractive theoretically. Some 
literatures have shown that the Kubo formulae with ladder resummation are 
equivalent to the kinetics theory in scalar field and gauge field [18,27,29,30]. 
In principle, one should resum various diagrams with ladder vertices and 
effective propagators when evaluating the correlation functions in Kubo for- 
mulae, which leads to a set of complicated integral equations. In this paper, 
we develop another approach, distinguished from the ladder resummation 
and the naive one-loop calculation which gives an inconvincible result, try- 
ing to demonstrate that even at one-loop skeleton diagram level the Kubo 
formula is still consistent qualitatively with kinetics theory if only the proper 
effective photon and fermion width are chosen as the infrared regulators in 
which the main physical characters are involved. 

In heavy-ion physics, the chemical potential in the central fire ball is not 
zero thought small compared to the extremely high temperature environ- 
ment. And most of the existing publications concentrated only on the high 
temperature case but neglected the small chemical potential in the heavy ion 
collision. To involve the chemical potential effect is our main motivation. 
It makes sense to introduce a small chemical potential <S T) to study 
that to what extent it can affect the properties of the plasma. The trans- 
port coefficients of QGP at finite density has been discussed in Ref. [22] with 
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the relaxation-time approximation, and early estimated by Danielcwicz and 
Gyulassy to scale the hydrodynamical equations [10]. We estimated the shear 
viscosity of hot QED plasma at finite temperature and chemical potential, 
referring to the results in the kinetics framework, since the infrared pinch 
singularities will appear and long distant interaction should be considered in 
the transport processes. 

The whole context will be arranged as follows: In section 2, we review 
the Kubo formula relevant to shear viscosity in the finite temperature theory. 
The contributions from the photons and the electrons in the QED plasma 
are studied respectively in section 3 and 4, in which one can see the trans- 
port damping rates play important roles in regulating the pinch singularities. 
Conclusions and discussions are presented in section 5. 

Here are some notations in this paper: (a) the capital letter stands for 
the four-momentum, K = (/cq, k) with k = |k|. (b) J is abbreviated to 
J dK. (c) When we mention 'hard', it means the momentum is much larger 
than T, or at least at the same order; 'soft' means much smaller than eT or 
at least at the same order. 



2 Shear Viscosity and Two-point Retarded 
Green Function 

Generally speaking, the longitudinal expansion domains in heavy ion colli- 
sion, therefore the shear viscosity is more important than the bulk viscosity 
which is proved much smaller than the former in lattice calculation [17]. In 
this paper, we ignore the bulk viscosity and just take the shear viscous one 
into account. 

In a near equilibrium system with linear response theory, Kubo formula 
tells us the shear viscosity [30] is determined by 

^ = ^j d'^' f dt'GRiO;x',t'), (1) 

where Gji{0;x',t') is the two-point retarded Green function defined as: 

Gr{x, t; x', t') = -id{t - t') < [TTy (x, t),nji{x', t')] >, (2) 
and TTij is the spacial part of the dissipative energy-momentum tensor 



TTjj = {SikSji - ^Sij6ki)T'''' = -Tj (^diUj - djUi + '^dijdku''^ , (3) 
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where T is the spacial energy momentum. The coefficients rj is the shear 
viscosity. The four-velocity u^{x) in the local rest frame is (1,0,0,0). 
Transforming Eq.Q into momentum space, one obtains 

with Gr{Q) = JdX'e-'^-^'GR{0;t',x') and X' = {t',x'). 

Through the Kubo formula we know that the shear viscosity is determined 
by the two-point retarded Green function. Now we start to evaluate it in an 
explicit dynamics. In the relativistic heavy ion collision, people are interested 
in the extremely high temperature environment, therefore one can use hard 
thermal loop(HTL) approximation theoretically. We here choose QED as a 
sample model and calculate the contributions both from the bosons and from 
the fermions. 

The QED Lagrangian reads: 

Cqed = ^{tl^ ■ d, - m)^ - ^F'^'F^, - e^Y'^A^, (5) 

where F^,^ = d^A,y — d^A^ and \E', A^ are the electron, anti-electron and 
photon fields respectively. 

With this Lagrangian, one can read the viscous stress tensors in the mo- 
mentum space, 

<^.(fc) = z^(7.fc,-l%7-k)^, (6) 

for electrons and 

i.r{k) = -{hk,-\5,,e)A' (7) 

for gauge boson with the Lorentz gauge condition. 

Combining Eqs.©, Q, ® and (j?!), with the standard definitions of 
propagator in quantum field theory and transport vertex defined behind, one 
can figure out the one-loop skeleton diagrams of the retarded Green function 
presented in Fig.l. The main difference between the propagators of the 
normal field and the thermal field rest on the fact that the latter is a 2 x 2 
matrix due to the doubled Hilbert space [31-33]. 

In a high-temperature environment, the loop momenta are hard therefore 
the mass of fermions can be omitted and bare thermal propagators are suf- 
ficient. Yet in the straightforward calculation pinch singularities occur when 
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(b) 



Figure 1: Two point retarded Green function of boson(a) and fermion(b). 
The dashed hne represents the spacial stress tensor operator TTjj, the sohd 
hne for fermions and the wiggle hne for bosons. 

the external four-momentum approaches to zero, which is required by Eq.(0}. 
To avoid such divergences, we must keep the imaginary part of the propaga- 
tor, namely, the effective damping rate(widths), as infrared regulators, which 
can be obtained analytically at leading-log order. 

3 Boson Contributions 

The diagram of the boson two-point Green function is demonstrated in 
Fig. 1(a), with which one can expressed it as 

GliiQ) = ^ I dKD^J^{K)Ilr,D,,^,,{K + g)/7r,, (8) 

where D^{Q) is the thermal gauge boson propagator, and a, b=l, 2 are the 
doubled Hilbert space indices; r^, = ±1 represent the type-1 and type- 
2 vertices respectively. The notation I]- = —kikj + ^Sijk"^ is the transport 
vertex of the photon viscosity. Summing over the Lorentz indices and using 
the reverse relation of Keldysh representation [31,34] we obtain: 

Gl\K) = GfAQ) - G',^{Q) (9) 
322 r 

= — I dKk^[nBiko) - nsiko + qo)]DA{K)DR{K + Q), 

in which nsiko) = (e^''^"' is the Bose-Einstein distribution. Dji^A{K) is 

the retarded(advanced) gauge boson propagator in the Feynman gauge with 
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D'S'JK) = g>"-DaAK) and 

In the last line of Eq.(jHI), the fluctuation-dissipation theorem has been em- 
ployed: 

Gs{P) = [1 + 2nB{po)]sgnipo)[DRiP) - Da{P)]. (11) 

Here we encounter the pinch problem of Da{K)Dji{K+Q) in Eq.Q when 
the external momentum Q approaches to zero which is required by the defi- 
nition of viscosity. The physical reason responsible for this divergence is the 
naive adoption of bare propagators. Therefore the so-called HTL resummed 
propagator developed by Pisarski et al. needs to be considered to involve 
the medium effect, namely, the scheme of eliminating the pinch is to endow 
the finite particle life-time or width consistently. Mathematically, we insert a 
small 7p/i into the imaginary part of the propagator which physically means 
the damping rate. The real part of self-energy(thermal mass term), which 
is in gT order, can be neglected compared with the hard loop momentum. 
We drop out Q in the propagators directly since Eq.(jll) requires go = after 
derivation on qq, and thus the thermal distribution in Eq.® has the form 

riBih) - nsih + go) ~ go/?ni3(A;o) [1 +^^(^0)] H (12) 

which preserves the first order in go. This procedure can be written manifestly 

as 

0.(K)D.(K,Q) (13) 

-^j^\S(k„-k)+S(k„ + k))]. (14) 

Notice that the photon damping rate is of order e^T , much smaller than the 
momentum K >T, which implies the width of the photon can be substituted 
by a delta function approximately in expression ()14|1 . Integrating off the delta 
function and the azimuth angles in the spherical coordinates, the retarded 
Green function becomes 

GtiQ) = ^ / dkk'nBim + Mk)] — . (15) 
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Figure 2: The physical processes contribute to the photon damping rate. 
The first fine is for the Compton scattering and the last line for the pair 
production. 



Now we are in the position to calculate the photon damping rate consis- 
tently at a finite temperature and density. There is a fact one should notice 
that the longitudinal contribution is suppressed compared to the transverse 
one for hard momenta [33]. To the leading order evaluation, one should only 
consider the latter. It was calculated by Thoma with vanishing chemical 
potential [35] , and we will generalize it to the finite density case. 

The photon damping in the plasma owes to two processes, Compton scat- 
tering and pair production(Fig.2). The two diagrams in the first line of fig. 2 
stand for Compton scattering and the last line is for pair production pro- 
cesses. 

The Braaten-Pisarski method [36] has been used to calculate the photon 

damping rate in the QED plasma, where the total effect is divided into two 
parts by a separating scale. The soft contribution is obtained by cutting 
the photon self-energy diagram with the HTL resummed fermion propaga- 
tor(Fig.3) and the hard contribution can be evaluated most conveniently from 
the matrix elements. Adding up the two contributions one can cancel the 
separating scale. 

3.1 Soft Contributions 

The soft contribution of a real hard photon damping rate follows from the 
the imaginary part of transverse self-energy in Fig.3, where the blob denotes 
the full fermion propagator. 
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Figure 3: Hard photon self-energy with the full propagators denoted by a 
blob. 



The damping rate of a photon is defined as 

7 = -^ImUripo = P,p), (16) 

in the case of non-overdamping. The self-energy is given by 

n^, = 2ie^ J dKTr[S*{K)^^S(P - K)^,] (17) 

where S*{K) is the effective fermion propagator, the factor 2 accounts for 
both diagrams in Fig. 3. We take Q = P — K and adopt the helicity repre- 
sentation of the fermion propagator, 

- 1 7o-k-7 I 1 7o + k-7 .-.ox 
^^^^ - D^{K) 2 ^D_{K) 2 ^^^^ 

^(O) - ^ 70 -q- 7 1 7o + q-7 (.a) 

where k = k/A; and 

D^{K) = -^o±fc + ^[(lT^)ln^±2] (20) 

d±iQ) = -go ±9, (21) 

with z = ko/k. Here the full fermion propagator is the result from HTL 
resummation. Projecting the self-energy on the transverse direction, one 
finds 

nT(Po,p) = ^(5,,-^)H,,(po,p). (22) 
z p 
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Evaluating the trace over 7 matrices, one can transform the transverse part 
of photon self-energy [37] into 

. . _ 2 /" ,^ /^l^ l + V \ 1 f l + V 1-V \ 

iVr{n - y dK ^^^^^^ y^^^^^ + J + y^^^^^ + J 

(23) 

where V = (p ■ k)(p ■ q). 

To sum the frequencies over the fourth momenta with a chemical po- 
tential, we first prove an auxilary formula which is useful to sum over the 
Matusubara frequencies in a typical one loop diagram at finite chemical po- 
tential. The zero chemical potential case was first developed by Braaten, 
Pisarski and Yuan [38]. 

The typical one loop diagram one encounters in the temperature field 
theory may take the form of 

0{iuj) = T A(iuJn + /i) A(ico' — itOn — /i) (24) 

n 

where iu is the external line momentum and iun = 27r(n + l)/f3 are the 
fermion Matusubara frequencies. Employing the spectral representation 



A{zuj^ + li) = - I dko . , , (25) 



one recasts Eq. ()24j) into 

0{iuj) = I dko [ dqop{ko)p'{qo)Ty^ \- 

J J ^ lUJn-k 



iun — ko + fi iu — iun — ko + fl 

(26) 

Summing over the Matusubara frequencies n, one obtains 

Oizuj) = - [dko [ dqop{ko)p'{qo)^—^^^^:^^^^ (27) 
J J lu-ko-qo 

with f±{ko) = [e^('=«T'^) + 

Let iuo be analytically continued to po which are not the Matusubara 
frequencies and the imaginary part of 0{po) is 

ImO(po) = ^BiscO{po) = ^^[0{po + te)-0{po-ie)] (28) 
= 7r(l-e'^P«) j dko j dqof+{ko)f-{qo)p{ko)p'{qo)5{po - ko - qo)). 
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When T » » /i, Uiko)f-{qo) = U{ko)f-iPo - h) ~ e-^^'V2, 

ImO(j9o) = |(1 - e^^'«)e-^^'» j dk, j dq,p{k,)p{q,)5{p, - k, - q,). (29) 

Using the generahzed Eq.(|^^ and adopting the similar steps in Ref. [37], 
we finally obtain 

^ 4p V 7r«(T2 + g) ^ ^ 

where a = e'^/Air is the fine structure constant, and A is the separation scale, 
namely, the upper limit of three-momentum integration. 

3.2 Hard Contributions 

The hard contribution, when the exchanged fermion carries momentum much 
larger than the separating scale A, can be evaluated through the matrix 
elements. The damping rate relevant to the Compton scattering process is 

x{27r)^6\P + K~P'-K')A<\M\'>,,mp, (31) 

where < \M\ >comp is the scattering amplitude and npik) = (e^'^"' + 1)^"*^ is 
the Fermi-Dirac distribution. Similarly, the damping rate for pair creation 
process is 

X {271)^ 6\P + K-P' - K')2 <\M\^ >pa^r ■ (32) 

Using Mandelstam variables s = {P + K)"^, t = {P — P'Y and u = —s — t, 
one finds the matrix elements are given by [39] 

<\M\'->comp = -2e' (- + -), (33) 

<\M\'>pair = 2e'(^+^-y (34) 
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In the hard external hnes assumption, that is P' ^ T and K' ^ T, 
the phase space for the outside region is unfavorable, thus the distribution 
functions can be simplified as 1 -^np^B ~ 1- This considerable simplification 
may bring us much convenience when using the Lorentz invariant phase space 
factor [39] 

dL = {27t)^S\P + K - P' - K" 



(27r)32p' (27r)32fc' 
* (35) 
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In our approximation, the density effect, namely the chemical potential, 
modifies only the Compton scattering process due to the initial fermion dis- 
tribution, while those in pair production are kept unchanged as in the pure 
temperature environment. And the result is already given by Ref. [35] 

7C^. = ^^(ln^- 2.1472). (36) 

In the Compton process, considering the initial fermion distribution function 
with a small chemical potential, we can obtain the corresponding damping 
rate 

where 

A{fi) = S2ie~^n + S2ie^n, (38) 

oo ^ 

'^«W = 5Z|^1^^' ri = 0,l,2---. (39) 



k=l 



Adding up Eqs. ljHBj) and (|H7|). one obtains the hard contribution. 



- 1. 2 



(T2 + ^) In ^ - (1.4572 + '^^)T^ - 0.0078/.^ 
A"' vr^ 



(40) 



When the soft and the hard contributions are combined together, the 
separation scale A is to be cancelled thus yielding 

T Isoft ~l~ Ihard 

= + (41) 
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with 

{p)T 



F(T,^) = ln. 



(lA572 + ^^^^ + 0.0078/1^ 



(42) 

In the above expression of F{T, /i) we replaced the momentum in the log- 
arithm by its average value (p) to simplify the future integration, which is 
defined as 

, , ( (Pk k nuik) 

Inserting the Eq. (j41|) into (jlSp . completing the three momentum integra- 
tion and with Eq.(j3]), we finally obtain: 

1 072T^ 

Vvh = —2 ~ 84.64— -(1 - 0.101^). (44) 

where the last approximation is obtained with expansion in terms of fi'^/T'^ 
and keeping solely the leading-log accuracy. 



4 Fermion Contribution 

In this section we are going to calculate the shear viscosity contributed from 
the electrons. With Eq. (jH)) and Fig. 1(b), one can write 

G:,iQ) = I j dKTr[Satim,nSatm (H )l]{fal (45) 

where If^ = 'jikj — |5ij7 ■ k and Sabi^P) is the element of 2 x 2 fermion 
propagator matrix. Summed over the thermal indices, the retarded Green 
function then becomes: 

Gl{Q) = Gn(g)-Gi2(Q) (46) 

Using the inverse relation of Keldysh representation and the fluctuation- 
dissipation theorem, the last line of Eq. P6|l recasts into 

z j dKTr[f{m^+(^)I]MMko) - npiko + qo)]Sn{K)SA{K + Q), (47) 
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K(E,k) 




Figure 4: Hard fermion damping in the QED plasma. The blob on the photon 
line stands by a resumed propagator. 



where hpiko) = [e^(l'=ol±'') + 1]-^ and S{f{) = fCS{K). 
We notice that 



and 



Sr{K)Sa{K + Q) ^ 



+ iesgniko) K"^ — iesgn{kQ 



(48) 



(49) 



Similarly, considering the fermion damping in the plasma and using the elec- 
tron damping rate 7e to improve the infrared behavior, one obtains 



Sr{K)Sa{K + Q) 



1 



1 



{ko + i7e)2 - k^ {ko - I'-feY - k 



2 ■ 



(50) 



The electron damping rate, defined by the imaginary part of the elec- 
tron self-energy on mass shell with full photon propagator (Fig.4), is just one 
half of the interaction rate T. According to the cutting rules [40,41], the 
interaction rate can be described by the photon longitudinal and transverse 
spectrum functions pi and pt defined in the HTL approximation for space-like 
frequencies as 



V{E) 



27rv 



POO P 

Jo 



vq ^ 



ixl 



pl{uJ,q) + (1 :r)pt{uj,q) 



(51) 
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with V = k/E and 




where = |e^T^(l + is the photon thermal mass. Performing inte- 

grals oidq and (ia; with quasi-static approximation, i.e. v -C 1, one finds that 
the longitudinal part is screened by the thermal mass while the transverse 
part suffers a logarithmic infrared divergency due to the absence of magnetic 
screening. Introducing an infrared cutoff e by hand, one can count the power 
dependency of the coupling constant 

V{E) ~ e^r(ln e-^ + constant) , (54) 

which makes the viscosity r) ~ e~^, namely, the reciprocal of quadratic cou- 
phng constant, but not quartic as kinetic theory suggusts [8,13,42]. This 
disagreement hints us to discover the defects in our previous calculation. 

Returning back to Fig.(l), one would find that only involving effective 
propagator is not enough, which means the rungs between the two propa- 
gators have been omitted. To include such diagrams we must sum up all 
these ladder diagrams. While that is a rather complex procedure and a set 
of integral equations is obtained as a result [18,27,29,30]. Here we do not 
employ this scheme, instead we replace the fermion damping rate, namely 
the infrared regulator, with a more physical transport damping rate, which 
has been applied even in the discussion of classic system and Abelian or 
non-Abehan plasma transport coefficients [10,21,22,24,26]. 

For a Rutherford-like scattering in classic plasma, namely, the elastic scat- 
tering with the same type and momentum of incoming and outgoing particles, 
small angle scattering domains in the dynamical cross section. Whereas, as 
we know, the large angle scattering, which is caused by a series of accumu- 
lated multiple small angle scattering in a long distance, is the most effective 
to the shear viscous process. Thus an additional factor of 1 — cos d is multi- 
plied in front of the ordinary cross section [43], where is the scattering angle. 
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In this procedure the so-called transport cross section is defined. Further- 
more, when discussing the transport process in QED or QCD plasma with the 
Boltzmann equation in kinetics theory, one will find that the collision term 
is the convolution integral of scattering amplitudes and their correspondent 
statistical weights, where if the interaction vertex connects the two particles 
of the same type, the fiuctuation part of the statistic distribution functions 
will contribute an extra small [13], 

[xij(p + q) - x^i(p)]' = [q ■ Vxi,(p) + ■■■]' (55) 

which is equivalent to multiply the cross section by the factor of 1 — cos 9, 
since 

l_cos^= — (1- i^). (56) 

s 

Weldon investigated the decay rate in thermal field theory [40], pointing 
out that the statistical distributions involving both direct and inverse reac- 
tions(detailed balance) should be considered simultaneously, which naturally 
demanded for introducing the transport cross section when calculating the 
fermion damping. 

Now we are aware of what has been lost in the straightforward compu- 
tation: that is q^ in Eq.(|3^. i.e., the statistical weights do not refiect the 
detailed balance, which leads to an incorrect result. Therefore the right way 
to obtain the effective fermion regulator (width) in transport process is to 
employ the transport interaction rate and yields, 

1 If 

le = -^trans = 7T / ^^(1 - COS 6^) (57) 



^2 

pi{uj,q) + (1 -)pt{uj,q) 



2 2„ 
and 

g2 roo rvq Jl 



trans 



2 r-oa p 

{E) = — dq q' / 

TTVS Jo J_ 



vq ^ (f 



(58) 

Changing integration variables from to x = tu/g, we can evaluate the 
integral over dq analytically and find out the logarithmic dependence on the 
integral limits. This is again the category which can be treated by employing 
the Braaten-Pisarski method for the consistent leading order quantities, i.e. 
employing a separating scale A to divide the integral into a soft contribution 
and a hard contribution, then adding the two contributions which leads to 
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the cancellation of A. One would notice in the previous calculation on the 
photon damping rate that the coefficients in front of the log-terms in the soft 
and hard contributions are identical so as to guarantee the separating scale 
to be cancelled. Thereby we can play a little trick to obtain the leading-log 
result by just calculating the soft contributions and replacing the separating 
scale with the upper T, ignoring the constant factor added to the logarithm 
term which is less important in the weak coupling limit. Completing the 
integral over dq we obtain^ 



7e 



2 ^ trans (-^) 



m log . 

Atis ' 



(59) 



with s = (P + KY = {'^EpY = in central mass system. 

Combining (jKlHl with (|77jl and noticing that 7e ~ e^T <^ T, one adopts a 
delta function as approximation and finds: 

GUQ) = / dK-[Mko) - Mko + qoWiko - k) + 6{ko + k))]. (60) 

3 i 7e 

Integrating over dko and then expanding the fii? in terms of small go, ignoring 
the terms without go (because Eq. (jH) needs go = after differentiation on 
go), we obtain 



GUQ) 



Svre^m^ In e ^ 



dkk^ 



(61) 



Noticing T ^ /i, one can expand the integrand of (jHH) as Taylor series 
with respect to fi/T and complete the final integration. The leading- log shear 
viscous coefficient contributed by the fermions is: 



r^e = 361.4- 



T3 



:Mni 



0.1765^). 



(62) 



^ Since the plasma is in a high temperature but small chemical potential environment, 
i.e. T ^ fj,, the fermion lines in figure 4 should be hard and the Pauli blocking has tiny 
effect on the final result. Compared with the contribution from Bose-Einstein distribution 
function of photon, the chemical potential in the Fermi-Dirac distribution can be neglected. 
As to the density effect in the full photon propagator, the small chemical potential just 
modifies the constant in the logarithm but preserves the order in aT. Qualitatively one 
can adopt this result. 
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On the contrary, if an interaction vertex connects the two particles of dif- 
ferent types, Uke Compton scattering or pair production we have discussed in 
the previous section, no kinetics factor will be isolated from the statistical 
weight, i.e., 1— cos 6 will not appear in the photon width expression [21,22,25]. 
That is why we got the right dependence of the coupling constant by just 
adopting the ordinary photon damping rate as the infrared regulator in the 
photon contribution. 

5 Conclusions and Discussion 

Adding up the contributions from both the electron and photon sectors, the 
shear viscous coefficient of QED plasma is obtained, 

riQED ^Ve + Vph = 446.0-^:^(1 + 0.1234^). (63) 

One may find that the density-relevant part of viscosity coming from the 
photons is negative but the one from electrons is positive. They compete 
with each other and leave a positive contribution of chemical potential, which 
demonstrates an enhancement effect of density. 

The above result is consistent, in the order accuracy of coupling constant, 
but differing by a factor with those obtained in kinetic theory at the finite 
density [42] and the pure temperature case [13] when //, goes to zero. What 
we mean consistency here is just at the qualitative level. More accurate effect 
like ladder resummation, in which none damping rate should be replaced by 
the transport one since it is self-consistently including the main characters of 
the transport processes, should be involved in the calculation. The one-loop 
approximation may partly account for the discrepancy of the factors scaled 
by the logarithms. 

Prom the proceeding analysis we can conclude that the estimation of the 
shear viscosity of QED plasma calculated at one-loop skeleton diagram level 
by choosing proper fermion width in the framework of the finite temperature 
field theory, is economical and reliable. Other transport coefficients of the 
relativistic plasma at finite temperature and density could be discussed simi- 
larly. Using this method helps us to handle the pinch singularity easily if the 
transport hard fermion and boson damping rates are known. Nevertheless, 
we may point out that some further improvements could be introduced to 
obtain more precise results such as the Landau-Pomeranchuk-Migdal(LPM) 
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effects due to tlie interference of multiple scattering process [44,45] for the 
complete leading-order contribution. 
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